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Abstract: We consider the Nambu-Goto bosonic string model as a description of the 
physics of interfaces. By using the standard covariant quantization of the bosonic string, 
we derive an exact expression for the partition function in dependence of the geometry of 
the interface. Our expression, obtained by operatorial methods, resums the loop expansion 
of the NG model in the "physical gauge" computed perturbatively by functional integral 
methods in the literature. Recently, very accurate Monte Carlo data for the interface free 
energy in the 3d Ising model became available. Our proposed expression compares very 
well to the data for values of the area sufficiently large in terms of the inverse string tension. 
This pattern is expected on theoretical grounds and agrees with previous analyses of other 
observables in the Ising model. 
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1. Introduction 

The idea that string theory may provide the effective description of confining gauge theories 
in their strong-coupHng regime is an old and well motivated one [1, 2, 3]; in this context, 
the string degrees of freedom describe the fluctuations of the colour flux tube. 

In the last years much effort has been devoted to test this conjecture. In particular, 
several results have been obtained in an "effective string"' approach, in which the conformal 
anomaly due to the fact that the theory is quantized in a non-critical value d 7^ 26 of the 
space-time dimensionality is neglected. In principle, this is a very problematic simplifica- 
tion, since conformal invariance is at the very heart of the quantization procedure (see the 
Conclusions section for some remarks on this issue). However, it was observed very early 
[4] that the coefficient of the conformal anomaly vanishes for large distances, i.e., for world- 
sheets of large size in target space. In fact, in recent years, thanks to various improvements 
in lattice simulations [5, 6, 7, 8] the effective string picture has been tested with a very 
high precision and confidence [8]- [28] by considering observables such as Wilson loops and 
Polyakov loop correlators. It turnes out that at large inter-quark distances and low tem- 
peratures the effective string, and in particular the simplest model, the Nambu-Goto one, 
correctly describes the Monte Carlo data. As distances are decreased, clear deviations from 
this picture are observed [11, 13, 17, 18]. 
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Within the effective string framework, one can fix a "physical gauge" (see section 1.2 
for some more details) and re-express a generic effective string model as a 2d (interacting) 
conformal field theory of d— 2 bosons. In this set-up, the inverse of the product of the string 
tension a times the minimal area A of the world sheet spanned by the string represents the 
parameter of a loop expansion around the classical solution for the inter-quark potential. 
The first term of this expansion yields the well known Liischer correction [2] . The second 
term was evaluated more than twenty years ago in [29], with a remarkable theoretical effort, 
for several different classes of effective string actions. Higher order corrections would require 
very complicated calculations and have never appeared in the literature. 

Recently, a set of simulations both in SU(2) and SU(3) lattice gauge theories (LGT's) 
[9, 10, 12, 17, 18, 20, 21, 22, 24] and in the 3d gauge Ising model [6, 8, 11, 13, 14], thanks to 
powerful new algorithms, estimated the inter-quark potential with precision high enough 
to distinguish among different effective string actions and to observe the contribution of 
higher string modes. To compare the effective string predictions with these new data in a 
meaningful way it is mandatory to go beyond the perturbative expansion. 

So far this has been done only for the simplest effective string action, the Nambu-Goto 
one [30], and for the cylindric geometry, which physically corresponds to the expectation 
value of the correlator of two Polyakov loops. In this case it has been possible to build 
the partition function corresponding to the spectrum of the Nambu-Goto string with the 
appropriate boundary conditions derived long ago in [31, 32], and to make a successful! 
comparison with the simulations [13, 14]. In [33] this partition function was re-derived via 
standard covariant quantization, showing that it indeed represents the exact operatorial 
result which re-sums the loop expansion of the model in the physical gauge. 

As a further step in this direction we derive in the present paper the exact partition 
function of the Nambu-Goto effective string for a toroidal world-sheet geometry, and we 
compare this prediction to a set of recently published [34] high precision Monte Carlo 
results for the corresponding observable in the 3d gauge Ising model, namely the interface 
expectation value. 

Indeed, the toroidal geometry corresponds in LGT's to the "maximal 't Hooft loop" 
(see, for instance, [15, 16]). However a much simpler, yet physically very interesting, 
observable can be associated to the same geometry if we consider a three-dimensional LGT 
with a discrete abelian gauge group (like the 3d gauge Ising model). In this case the 
gauge model is mapped by duality into a three-dimensional spin model (like the 3d spin 
Ising model). In particular, the confining regime of the gauge model is mapped into the 
broken symmetry phase of the spin model. Any extended gauge observable, like the Wilson 
loop or the correlator of two Polyakov loops, is mapped into a set of suitably chosen anti- 
ferromagnetic bonds. The torus geometry we are interested in corresponds to the case in 
which the set of anti- ferromagnetic bonds pierces a complete slice of the lattice, i.e., to 
the case in which we simply impose anti-periodic boundary conditions in one of the lattice 
directions. In a spin model with discrete symmetry group this type of boundary conditions 
is known to create, in the broken symmetry phase, an interface between two different vacua. 
The (logarithm of the) expectation value of the slice of anti-ferromagnetic bonds, which we 
expect to be described by the Nambu-Goto string, is thus proportional to the interface free 
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energy, an observable which has been the subject of several numerical and experimental 
studies in condensed matter literature. We shall briefly recall these results in section 1.1. 

To compute the partition function, we follow the philosophy of [33] and resort to 
standard covariant quantization in the first order formulation of the Nambu-Goto theory, 
of which we briefly recall some aspects in section 1.2. In section 2 we integrate appropriate 
sectors of the bosonic string partition function over the world-sheet modular parameter. In 
this way we describe in an exact manner the string fluctuations around a specified target 
space surface, in our case representing an interface in a compact space. In section 3 we 
show that our expression reproduces the result obtained in [29] via a perturbative expansion 
(up to two loops) of the NG functional integral; in fact, our expression re-sums the loop 
expansion. 

In section 4 we compare our predictions with the data for the free energy of interfaces 
in the 3d Ising presented in [34]. Our NG prediction agrees remarkably to the data for 
values of the area larger than (approximately) four times the inverse string tension, which is 
the same distance scale below which deviations from the NG model emerged in the studies 
of other observables cited above [13]. In the range where there is agreement, we find that 
the NG result is largely dominated by the lowest level mode of the bosonic string; the free 
energy associated to this mode already accounts very well for the MC data, apart from a 
shift of the overall normalization. Though corresponding to a single particle mode, this 
contribution is essentially stringy, since this mode pertains to a wrapped string. 

1.1 Interfaces and effective strings 

The properties of interfaces in three-dimensional statistical systems have been a long- 
standing subject of research. In particular the interest of people working in the subject has 
been attracted by the so called "fluid" interfaces whose dynamics is dominated by massless 
excitations (for a review see for instance [35, 36]). For this class of interfaces, thanks to 
the presence of long range massless modes, microscopic details such as the lattice structure 
of the spin model or the chemical composition of the components of the binary mixture 
become irrelevant and the physics can be rather accurately described by field theoretic 
methods. 

An effective model widely used to describe a rough interface is the capillary wave 
model (CWM) [37, 38]. Actually this model (which was proposed well before the Nambu- 
Goto papers) exactly coincides [36] with the Nambu-Goto one, since it assumes an effective 
Hamiltonian proportional to the variation of the area of the surface with respect to the 
classical solution. 

A simple realization of fluid interfaces is represented by 3d spin models. In the broken- 
symmetry phase at low temperature, these models admit different vacua which, for a suit- 
able choice of the boundary conditions, can occupy macroscopic regions and are separated 
by domain walls which behave as interfaces. For temperatures between the roughening and 
the critical one, interfaces are dominated by long wavelength fluctuations (i.e. they exactly 
behave as fluid interfaces); all the simulations which we shall discuss below were performed 
in this region. 
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In these last years the 3d Ising model has played a prominent role among the various 
realizations of fluid interfaces, for several reasons. The universality class of the Ising model 
includes many physical systems, ranging from binary mixtures to amphiphilic membranes. 
Its universality class is also the same of the 0^ theory; this allows a QFT approach to the 
description of the interface physics [39, 40, 41]. Last but not least, the Ising model, due to 
its intrinsic simplicity, allows fast and high statistics Monte Carlo simulations, so that very 
precise comparisons can be made between theoretical predictions and numerical results. 

Following this line, during the past years some high precision tests of the capillary 
wave model were performed [42, 43, 34]. A remarkable agreement was found between 
the numerical results for interfaces of large enough size and the next to leading order 
approximation of the CWM. The present paper represents a further step in this direction, 
since we shall be able to compare the numerical data of [34] with the exact prediction of 
the Nambu-Goto effective model for the interface free energy. 

1.2 The Nambu-Goto model 

Perhaps the most natural model to describe fluctuating surfaces is the Nambu-Goto bosonic 
string [30], in which the action is proportional, via the string tension o", to the induced 
area of a surface embedded in a d-dimensional target space: 

(f^^d^ , gap = -Q^-Q^Gij . (1.1) 

Here we parametrize the surface by proper coordinates and {i = 1, . . . ,d) de- 

scribes the target space position of a point specified by ^. For us the target space metric 
Gij will always be the flat one. 

The invariance under re-parametrizations of the action (1.1) can be used to fix a so- 
called "static" gauge where the proper coordinates are identified with two of the target 
space coordinates, say and X^. The quantum version of the NG theory can then be 
defined through the functional integration over the (i - 2 transverse d.o.f. X{X^,X^) of 
the gauge-fixed action. The partition function for a surface E with prescribed boundary 
conditions is given by 



(1.2) 



= y DX' exp I y dX^dX^ {\ + {d^Xf + {diXf + {doX A diXf^ ' 

= / DX' exp \-a [ dX^dX^ 1 + ^{doXf + ^idiXf + interactions I 
J(9S) I is [2 2 J J 

Expanding the square root as in the second line above, the classical area law exp(— c7„4) 
(where A is the area of the minimal surface S) is singled out. It multiplies the quantum 
fluctuations of the flelds X, which have a series of higher order (derivative) interactions. 
The functional integration can be performed perturbatively, the loop expansion parameter 
being l/{aA), and it depends on the boundary conditions imposed on the fields X, i.e., 
on the topology of the boundary and hence of S. The cases in which S is a disk, a 
cylinder or a torus are the ones relevant for an effective string description of, respectively, 
Wilson loops, Polyakov loop correlators and interfaces in a compact target space. 
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The computation was carried out up to two loops in [29], see also [42, 44]. The surface 
S is taken to be a rectangle, with the opposite sides in none, one or both directions being 
identified to get the disk, the cylinder or the torus topology. At each loop order, the result 
depends in a very non-trivial way on the geometry of S, namely on its area A and on the 
ratio u of its two sides: it typically involves non-trivial modular forms of the latter. The 
two loop result for the case of interfaces is reported here in section 3. 

An alternative treatment of the NG model, which is the standard one described in 
most textbooks on string theory, see for instance [45], takes advantage of the first order 
formulation, in which the action is simply 

S = aJdfJ^ d^^ h'^^daX'dpX' . (1.3) 

Here hap is an independent world-sheet metric, G [0, 27r] parametrizes the spatial exten- 
sion of the string and its proper time evolution. Integrating out h, we retrieve the NG 
action eq. (1.1). 

As it is well known, to include the process of splitting and joining of strings, i.e., 
to include string interactions, one must consider world-sheets of different topology, i.e., 
Riemann surfaces of different genus g. The genus of the world-sheet represents thus the 
loop order in a "string loop" expansion. For each fixed topology of the world-sheet, instead 
of integrating out h, we can use re-par ametrization and Weyl invariance to put the metric 
in a reference form e'^hajs (conformal gauge fixing). For instance, on the sphere, i.e. at 
genus (7 = 0, we can choose h^p = rya^, while on the torus, at genus 5 = 1, is constant, 
but still depends on a single complex parameter r, the modulus of the torus; see later for 
some more details. The scale factor e"^ decouples at the classical level^ and the action takes 
then the form 

//•27r 
j di^ k'^da.X'dpX' + 5gh. , (1.4) 

where S'gh. in eq. (1-3) is the action for the ghost and anti- ghost fields (traditionally 
called c and b) that arise from the Jacobian to fix the conformal gauge; we do not really 
need its explicit expression here, see [47] or [45] for reviews. The modes of the Virasoro 
constraints T^p = 0, which follow from the /i'*^ equations of motion, generate the residual 
conformal invariance of the model. The ghost system corresponds to a CFT of central 
charge Cgh. = —26. The fields X*(t, o"), with i = 1, . . . , d, describe the embedding of the 
string world-sheet in the target space and form the simple, well-known two-dimensional 
CFT of d free bosons. 



2. The partition function for the interface from bosonic strings 

In the present section we want to describe the fluctuations of an interface in a toroidal 

^Actually, this property persists at the quantum level only if the anomaly parametrized by the total 
central charge c = d — 26 vanishes; Weyl invariance is otherwise broken and the mode 0, as shown by 
Polyakov long ago [46], has to be thought of as a field with a Liouville-type action. However, as argued 
in the Introduction and in the Conclusions section, we can, in first instance, neglect this effect for our 
purposes. 
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target space T'^ by means of standard closed bosonic string theory. We use the standard 
first order formulation discussed in section 1.2 and we specify the periodicity of the target 
space coordinates to be ~ + {i = 1, . . . ,d). 

The partition function for the bosonic string on the target torus T*^ is expressed as 



= [ ^Z^''\q,q)Z^\q,q). (2.1) 

J T2 



Here r = ri + ir2 is the modular parameter of the world-sheet, which is a surface of genus 
g = 1, i.e., a torus. Moreover, Z^^\q,q) is the CFT partition function of the d compact 
bosons X* defined on such a world-sheet. In an operatorial formulation, this reads 

Z^'^\q,q) = Tr g^""^ g^O"^ , (2.2) 

where 

q = exp(27riT) , q = exp(— 27rir) . (2-3) 

Lq and Lq (particular modes of the Virasoro constraints) are the left and right-moving 
dilation generators. With Z^^{q,q) we denote the CFT partition function for the ghost 
system, defined on the same world-sheet. 

The modular parameter r is the Teichmuller parameter of the world-sheet surface: as 
discussed in subsection 1.2, using Weyl invariance and diffeomorphisms wc can choose the 
reference metric hap to be constant and of unit determinant, but the complex parameter r, 
with Imr > 0, which characterizes the complex structure, i.e., the shape of the torus, cannot 
be fixed. It is thus necessary to integrate over it, as indicated in eq. (2.1), the integration 
domain^ being the upper half-plane. In the Polyakov approach [46], this integration is the 
remnant of the functional integration over the independent world-sheet metric hajs after 
the invariances of the model have been used as described above. The measure used in eq. 
(2.4) ensures, as we will see, the modular invariance of the integrand. 

The CFT partition function for a single boson defined on a circle 

x{e,e)-x{e,e)+L (2.4) 

is given, with the action defined as in eq. (2.1), by^ 

Ziq,q) = TV g^O"^ g^^"^ = V qsbi'ir+-^Lf qs^{'-¥-^l^f -1- . (2.5) 

The integers n and w are zero- modes of the field X, describing respectively its discrete 
momentum p = 2TTn/L and its winding around the compact target space: X must be 
periodic in 4^, but the target space identification eq. (2.4) allows the possibility that 

xie,e+'27r)=x{e,e)+wL . (2.6) 



^We will discuss later the issue of discrete modular transformations acting on r. 

^With an abuse of notation, and for the sake of convenience, we will sometimes denote the Dedekind eta 
function ij{t) defined in eq. (C.l) as ri{q), where q = exp(27i-ir). 
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The factors of l/r]{q) and l/rj{q) result from the trace over the left and right moving non- 
zero modes, which after canonical quantization become just bosonic oscillators contributing 
to Lq and Lq their total occupation numbers. 

The Hamiltonian trace eq. (2.5) can be re-summed a la Poisson, see eq. (C.IO), over 
the integer n, after which it becomes 




-%\ra-Tw\^ 



(2.7) 



an expression which is naturally obtained from the path-integral formulation. The dis- 
crete sum over m,w represents the sum over "world-sheet instantons", namely classical 
solutions of the field X which, beside their wrapping number w over the direction, are 
characterized also by their wrapping m along the compact propagation direction: 



X{e + 2TTT2,e + 27rri) = Xie,^) + mL . 



(2.8) 



The form eq. (2.7) of Z{q,q) makes its modular invariance manifest. In fact, the 
combination y/T2r]{q)r]{q) is modular invariant, as it follows from the properties of the 
Dedekind eta function given in eq. (C.4). Moreover, from the exponential term we infer 
the effect of modular transformations of the parameter r on the wrapping integers w,m: 
they act as SL(2,Z) matrices on the vector {m,w). In particular, the S and T generators 
of the modular group are represented as follows: 



1 

r 

r + l , 




(2.9) 
(2.10) 



This allows to reabsorb the effect of modular transformations by relabelling the sums over 
m and w. 

The partition function for the ghost system is given by 



Z^\q,q) = {v{qMq)f , 



(2.11) 



namely it coincides with the inverse of the non zero-mode contributions of two bosons. 
Notice that the string partition function eq. (2.1) can be rewritten, substituting the above 
expression for Z^^{q,q), in an explicitly modular-invariant way: 



(2.12) 



so that the integration over the modular parameter r in eq. (2.12) has to be restricted 
to the fundamental cell of the modular group. Indeed, the Poincare measure d^r/(T2)^ 
and the combination \/T2r]{q)ri{q) are modular invariant. The bosonic partition function 
Z^*^) ((?, q) is also modular invariant, being the product of d expressions of the type eq. (2.5). 
It depends on integers and u;*, the discrete momentum and winding number for each 
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Figure 1: The mapping of the toroidal string world-sheet, of modular parameter r, into the target 
space is organized in many distinct sectors, labeled by the integers Wi and rrii (see the text). By 
selecting the sector with (say) wi = I and TO2 — 1 we are considering the fluctuations of an extended 
interface, which is a torus because of the target space periodicity. 
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Figure 2: Consider an embedding of the world-sheet torus into the target space aligned along 
the directions x^, a;^, characterized by the wrapping numbers (m^, w^) and (m^, w^). It corresponds, 
in the covering space of this T^, to a parallelogram defined by the vectors w — (w^Li,k;^L2) and 
m = {m^Li,m^L2), whose area is w Am = LiL2{'w^m^ — vP'tv}). 



direction. For each direction i, we can Poisson re-sum over the discrete momentum n*, as 
in eq. (2.7), replacing it with the topological number m*. 

We want to single out the contributions to the partition function eq. (2.4) which 
describe the fluctuations of an interface aligned along a two-cycle inside T'^, say the 
one in the x^,x^ directions. The world-sheet torus parametrized by must be mapped 

onto the target space torus by embedding functions X^(^'^, ^^), X^(^'^, ^^) with non-trivial 
wrapping numbers {m},w^) and {rn?',w^). The wrapping numbers of X'''{SP,^^), i > 2, 
must instead vanish. The minimal area spanned by such a wrapped torus is given by (see 
Fig. 2) 

m^n?) = L1L2 (ml wA \~}\ I "\ 1 • (2.13) 



L1L2 iw^vn^ 



1 




A particular sector which contributes to such a target-space configuration can be ob- 
tained (see Figure 1) by considering a string winding once in, say, the direction: 



W2 = W3 



Wd = 



(2.14) 
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Figure 3: On the left, (some of) the coverings corresponding to the chosen sector of the partition 
function are depicted. Beside the one which exactly corresponds to the fundamental cell of the 
target space T^, i.e. the case with (m^,w^) = (0,1) and (m^,w^) — (1,0) there are "slanted" ones 
corresponding to generic values of m^. The generators S and T of the world-sheet modular group, 
which act by SL(2, Z) matrices as indicated in eq. (2.9), map these coverings to different ones with 
the same area. For instance, on the right, we draw the S- and T-transform of the "fundamental" 
covering discussed above (the solid one in the leftmost drawing). 



and, upon Poisson re-summation along the directions to x"^, selecting the integers 

7712 = 1 , = = . . . = md = . (2-15) 

This corresponds to wrapping numbers in the x^, directions given by {m^,w^) = {m^, 1) 
and {m'^,w'^) = (1,0) and hence, according to eq. (2.13), to a minimal area L1L2: such 
a configuration covers the two-cycle in the x^,x^ directions just once. We have not fixed 
the value of the wrapping number 772-^, therefore the sector we consider contains infinite 
embeddings that corresponds to "slanted" coverings, see Fig. 3. 

The choice Wi = nii = iov i = 3, d is preserved under the action eq. (2.9) of 
the modular group. Moreover, this SL(2, Z) action preserves the area eq. (2.13), as it is 
easy to see^. The orbit of the modular group which contains the configurations chosen 
in eq.s (2.14-2.15) spans thus infinite other string configurations characterized by different 
wrapping numbers in the directions x^, x^; all these configurations, however, always wrap 

*The expression of the area eq. (2.13) contains the simplectic product of the two vectors {m^,w^) 
and (m^,w'^). The simultaneous action of SL(2,Z) on both vectors preserves this symplectic product: 
SL(2) ~ Sp(l). 
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the in target space just once, and correspond to equivalent descriptions of the interface; 
see Fig. 3. 

To obtain the partition function for this interface we should include all the wrapping 
numbers {m^,w^) and {m^jW^) which are in the modular orbit of the configurations in eq.s 
(2.14-2.15), but we should integrate the r parameter in eq. (2.16) over the fundamental cell 
of the modular group only, as usual, to avoid overcounting. We can, equivalently, consider 
only the configurations of eq.s (2.14-2.15), fixing in this way the degeneracy associated to 
modular transformations, and integrate r all over the upper half-plane, which contains all 
the images under modular transformations of the fundamental cell. 

The second choice turns out to be convenient, as the integrals become very simple. In 
this way, we get the following expression for the interface partition function: 



-I d-2 



Expanding in series the Dedekind's rj functions: 




Li e . 



(2.16) 



2-d 



E 

fe=0 



and expressing q and q in terms of r, we rewrite eq. (2.16) as 

/ / _ \ °° foo re 

= n (V^i.) E E <^'* / <ine-'('-''«.) / 

VV^TT J-oo Jo 



dro 



X exp < —T2 



aL 



fc,fc'=OnieZ 

27r2n2 



(ra) 2 



^- + —TT- + Mk + k' 
2 crLj 



d-2 
12 



1 






T2 


2 


} 



(2.17) 



(2.18) 



The integration over T2 can be carried out in terms of modified Bessel functions using the 
formula 



(2.19) 



with 

and 

£ -- 



A = 



Li £ 



47r ,, d-2, 47r2n? 
— {k + k' - -^^) + i 



1 + 



B 



(2.20) 



aL 



12 



a^Lf 



/ 47r ti , , , , d-2. 47r2 nl: 
1 + — r(fc + fc'--rT^) + 



aA 



12 



(2.21) 

In the second step, wc introduced the area A and the modular parameter iu of the torus 
swept out by the string in the target space, namely the interface: 



A = L1L2 , u = 



Li 



(2.22) 
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The integration over n produces a 6{k — k' + ni) factor, and thus, using Eq.s 2.19-2.20, 
we rewrite the interface partition function eq. (2.18) as 



X(<^) = 2(^) ' FtV^ J2 w(-) K^{aA£) , (2.23) 

^ k,k'=0 ^ 

where f , introduced in eq. (2.21), is now to be written replacing k — k' for m. In eq. 
(2.23) we have introduced the transverse volume Vr = 11^=3 -^i> have re-written the 

prefactor y/aL2 as \/aAu. 

Eq. (2.23) gives the exact expression for the fluctuations of the interface aligned 
along the x^,x'^ directions, in our hypothesis that they be described from standard bosonic 
string theory. This expression depends only on the target space geometric data, namely 
the transverse volume Vr and the area A and the shape parameter u of the interface. 

Considering the leading exponential behaviour of the Bessel functions only, eq. (2.23) 
would reduce to a partition function constructed by summing over closed string states, 
characterized by the left- and right-moving occupation numbers k and k' and by the mo- 
mentum rii which is fixed to k — k' by the level-matching condition, with the usual bosonic 
multiplicities c^, Cfc' and with exponential weights exp(— L2-E'A;,fe')) where the energies are 
given by 

/ rl —2 47r2 tj^ 

Ej,,k' = cjL^£ = JcT^Ll + 47ra(fe + k' - + . (2.24) 

This spectrum substantially agrees with the expression proposed in [48]; however, the 
correct expression of the partition function, eq. (2.23), involves Bessel functions rather 

than exponentials and it also contains extra factors of £. These modifications are crucial 
in reproducing correctly the loop expansion of the functional approach, as we shall see in 
the next section. 



3. Comparison with the functional integral approach 



As discussed in the Introduction, the Nambu-Goto interface partition function, whose 
expression eq. (4.3) we derived by opcratorial methods in the first-order formalism, can 
also be computed in a functional integral approach with a physical gauge-fixing. This leaves 
just the d — 2 bosonic d.o.f. corresponding to the transverse fluctuations of the interface, 
but the action is not a free one. One can evaluate the path integral perturbatively, the 
loop expansion parameter being the inverse of a A [29]. The result of this computation up 
to the 2-d loop order was given in [29] and reads 



:.-(tA 



where the dots stand for higher loop contributions and 



1 + ^ + ... 



(3.1) 



id - 2f 



+ 



d{d - 2) 



(3.2) 
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Here E2 denotes the 2nd Eisenstein series (see Appendix C). Actually, the constant term 
in eq. (3.2) is different from the one given in [29]; in Appendix B we show, by reconsidering 
the computation of the second loop term, that in fact the correct result is the one we quote 
ineq. (3.2)^ 

Our exact expression eq. (2.23) should reproduce the pcrturbative expansion of the 
functional integral result when asymptotically expanded for large a A. To this effect, let us 
re-organize eq. (2.23) in a suitable way. Using the asymptotic expansion of Bessel functions 
for large arguments: 

vr 



K,(z)~ J-e- l + ^a(^VM , (3.3) 



(7) 

where the coefficients of are well known, we obtain 



^ k,k'=0 \ r=\ J 

We can then Taylor expand the expression eq. (2.21) of £ and write it in the form 

CO 

S = l + ^dt{a,b-u){(jA)-\ (3.5) 
t=i 

having introduced, for notational simplicity, 

a = k + k'- , b = k-k' . (3.6) 

In particular one has 

di = 27rna = 2ttu (k + k' - ^^r^ ) . (3.7) 



12 

For the sake of readability, we leave to Appendix A most of the details of the present 
computation, such as explicit expressions of the various expansion coefficients, describing 
here only the basic steps. 

Plugging the expansion eq. (3.5) into eq. (3.4) we obtain 

x(^)oc.^^ f c.c,,e-M^+^'-^)|l + f;^i/^| , (3.8) 

where the coefficients Qs can be re-constructed starting from the coefficients appearing in 
eq.s (3.3) and (3.5); see Appendix A for more details. 

We can introduce Q = exp(— 27ru) and re-write eq. (3.8) as follows: 

mn f; J 1 + f; \ cf-^Q"'-'^ . (3.9) 



^It is interesting to observe that the missing contribution in [29] is proportional to (d — 3) and thus 
disappears in three dimensions. This is the reason for which it was not found in the calculations reported 
in [42, 44] which evaluated this 2 loop correction with three other types of regularization in the d — 3 case. 
Moreover the difference is modular invariant and thus it could not be detected by Dietz and Filk in the 
tests of their calculation which they made in [29]. 
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At this point, we can effectively replace 



, , d — 2 ^ d d 



b = k-k' —.Q^-Q^ (3.10) 
aQ dQ 



and take into account the definition, eq. (2.17), of the coefficients Ck to obtain the form 

-a A 



r,^. d-2 e 



U 2 Q^Q 



4^9s(QjQ-Qj^,Q^ + Q^;u)] _ 

1 + ^ -^-^ '^-^ \ [r){QMQ)] " 



s=l 



{aAy 



(3.11) 

with, in the end, Q being set to exp(— 27r'u). Finally, we can rewrite the above expression 
as 



-a A 



[^/urj'^{m)] 
where the coefficients fs are given by 



d-2 



fsju) 



(3.12) 



)_ hm g(Q±^-Ql^^Q±^+Q^.^[r^[Q)r^[Q)f-'' (3.13) 
Q-»Q=e-2''" \ dQ dQ dQ dQ J 



and can be related to Eisenstein series. In Appendix A we derive the expression of the 
2-nd loop and 3-rd loop coefficients, /i and /2, but it would be straightforward to push the 
computation to higher orders. Our expression for /i agrees completely with the (corrected) 
Dietz-Filk result cq. (3.2). Let us remark that our result eq.(3.12) is proportional (with no 

need of ad hoc modifications) to , i.e. to ^/a in three dimensions. This result agrees 

with the (completely independent) field theoretical calculations in the framework of the 
model of [39, 40] and also, as we shall see, with the Monte Carlo simulations. 



4. Compcirison with Monte Carlo data 

In a very recent publication [34] , precise Monte Carlo data on the free energy Fg of interfaces 
in the 3d Ising model were presented. The most accurate sets of data were obtained for 
square lattices (i.e., in the case u = 1), for values of the Ising coupling given by /3 = 0.226102 
(Table 4 of [34]; we shall refer to this data set as the set n. 1) and (3 = 0.236025 (Table 
2 of [34], set 2). The data set 1 contains 21 points, obtained at lattice sizes L\ — L2 = L 
ranging from L = 10 to L = 30; the 9 points in data set 2 are^ for L = 6 to L = 14. 

Previous work regarding other observables [13] has made it clear that, in appropriate 
regimes, the 3d Ising model can be successfully described by an effective string theory. 
For the two values of /? which we study here very precise estimates for the string tension 
exist [49] (see Table 1 of [34]). For set 1 one has a = 0.0105241; for set 2, cr = 0.044023. 
This entails that the points in data set 1 correspond to values of s/oA ranging from 1.02 
to 3.07. while (ho points in data set 2 correspond to values from 1.26 to 2.94. 

®We consider here only the data in Table 2 of [34] which, for a given value of L, are obtained with Lq, 
the size of the lattice in the transverse direction, being equal to 3L, to ensure uniformity with the data set 
1, which is obtained with this choice. 
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100 







-^min 


(\/cr^)min 


N 


xV(d.o.f) 


N 


xV(d.o.f) 


Data set 1 


19 


1.949 


0.91957(18) 


4.22 


0.91413(18) 


1.60 


20 


2.051 


0.91891(22) 


1.84 


0.91377(22) 


0.88 


21 


2.154 


0.91836(27) 


0.63 


0.91344(27) 


0.47 


22 


2.257 


0.91829(33) 


0.70 


0.91354(33) 


0.50 


23 


2.359 


0.91797(45) 


0.63 


0.91339(45) 


0.53 


24 


2.462 


0.91762(57) 


0.57 


0.91316(57) 


0.55 


25 


2.565 


0.91715(75) 


0.50 


0.91279(75) 


0.55 


Data set 2 


9 


1.888 


0.91052(21) 


7.22 


0.90466(21) 


2.22 


10 


2.098 


0.90924(33) 


2.71 


0.90413(33) 


1.69 


11 


2.308 


0.90820(51) 


1.12 


0.90349(51) 


1.33 



Table 1: The fit of the NG free energy eq. (4.1) with normahzation M to the two data set of Ref. 
[34], performed using only the points in Table 4 and 2 of [34] corresponding to lattice sizes L > Lmin, 
i.e.. those with \faA > (\/ov4)min- The reduced becomes of order unity for {\J oA)rmxi ^ 2. The 
fit is performed by truncating the sum over the oscillator levels at iV = 100 or at = 0, i.e. keeping 
only the 0-mode contributions (last two columns). 



Using the information above, it is possible to compare the MC values of the free energy 
Fs in data set 1 and 2 to the free energy F corresponding to our partition function eq. 
(2.23) (in d = 3, and for u = 1), where we factor out the transverse volume factor Vt and 
allow for an overall normalization e"-'^: 

The constant will be the only free parameter to be fitted to the data. For numerical 
evaluation of the free energy eq. (4.1), we re-arrange the sums appearing in eq. (2.23) 
using the property 

oo oo m 

""kCk' f{{k - k'f, k + k') = J2Yl ^kCm-kf {{2k - mf, m) , (4.2) 

k,k'=Q m=0 fe=0 

valid for any function / of the specified arguments, and write eq. (2.23) as 

d-2 OO m / p\ 

1^'^ = 2(1^) ' Vt J2 ^^^-^ ( - ) (^•^) ' (^-3) 

m=0 k=0 ^ ^ 

where £ is now £{{2k — m)^,m), since we transformed the summation indices as in eq. 
(4.2). The integer m represents the total (left- plus right-moving) oscillator number; the 
contributions to the partition function from states of a given m are suppressed by exponen- 
tials of the type exp(— 27r'um), as it can easily be seen from eq. (4.3) upon expanding the 
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Figure 4: Two sets of Monte Carlo data for the interface free energy data provided in [34] are 
compared to our tlicorctical predictions following from eq.s (4.1,4.3), represented by the solid red 
line for the data set 1 and by the dashed blue line for data set 2. The only free parameter is the 
additive constant A/", corresponding to an overall normalization of the NG partition function, fitted 
to the data using the points to the right of the vertical grey dashed line (see the text for more 
details). The error bars in the MC data are visible only for the rightmost points, but they were 
kept into account in the fit. The sum over the level m in eq. (4.3) was truncated at = 100. 

argument of the Bessel function using eq.s (3.5,3.7). We can approximate the exact expres- 
sion of F by truncating the sum over m at a certain N, i.e., neglecting the contributions 
of string states with total occupation number larger that N. 

Having fixed u = 1 and having specified A'^, the free energy F of eq. (4.1) is just a 
function of aA, which we can fit to the MC points in the two data sets referred to above. 

The results of this comparison are summarized in Table 1. For each data set, we 
consider only the data corresponding to lattice sizes L > Lmim i-e., to values of \/ a A > 
(V^) miri; and fix the value of the normalization constant AA so as to minimize the x^'i we 
repeat the analysis for various values of Lmin- Also, we consider the case where we truncate 
the sum over the oscillator number at = 100 and the case = where we keep only 
the zero-mode contributions. 

For both data sets 1 and 2, the reduced becomes of order one when we consider 
in the fit only values of V a A > 2. This is therefore the region in which our theoretical 
expression, derived from the Nambu-Goto model, describes well the data. This pattern 
is fully consistent with previous analysis of other observables in the 3d Ising models. In 
particular, the case of Polyakov loop correlators was recently discussed in [13], where 
two different behaviours were observed. In the compact direction corresponding to the 
"temperature" , a good agreement with the Nambu-Goto predictions was found exactly in 
the same range of distances as in the present case. On the contrary, in the direction with 
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Dirichlet boundary conditions a clear disagreement was observed up to very large distances. 
This is probably due to some type of "boundary effect" associated to the Dirichlet boundary 
conditions. It is exactly the absence of this type of effects which makes the case that 
we study in the present paper, periodic in both directions and hence free of boundary 
corrections, particularly well suited to study the effective string contributions, disentangled 
from other spurious effects. 

Considering the values of the we can take as our best estimates for the normalization 
constant J\f the value obtained at Lmin = 21 for the data set 1, namely J\f = 0.9184(3), 
and the value at Lmin = H for data set 2, i.e. M = 0.9082(5). It is interesting to observe 
that truncating to N = Q the mode expansion gives comparable values for the but 
leads to a sizeable change (with respect to the statistical errors) of which becomes 
N = 0.9134(3) and J\f = 0.9035(5) for the data sets 1 and 2 respectively. This difference 
is due to the fact that the truncation of the sum at A'^ = replaces the area-independent 
1-loop determinant ry'^^^''(iu) appearing in the exact partition function eq. (3.12) with its 
approximation exp(7r'u(d — 2)/6). This can be seen by keeping the m = 0,k = term 
only in eq. (4.3), looking at the exponential behaviour e^p{—aAS) and expanding S to 
the first order, see Eq.s (3.5-3.7). In our situation, where d = 3 and n = 1, the 1-loop 
determinant (which is essentially reproduced summing the contributions up to = 100) 
contributes to the free energy a constant term 21og[?7(i)] = —0.527344... ; in the A'^ = 
truncation, this constant becomes — 7r/6 = —0.523599. The difference between these two 
values is 0.00374536..., and accounts for most of the difference between the estimates of 
the overall constant J\f &t N = 100 and N = given above. 

It is instructive to compare our estimates for the overall constant with those obtained 
in [34] by truncating the perturbative expansion in 1/aA to the second order. Keeping 
into account the difference in the definition of the normalizations, which simply amounts 
to the contribution 2 log [?7^(i)] of the one-loop determinant, we would obtain, using the 
notations of [34], co + ^ log a = 0.3910(3) for the data set 1 and N = 100, which should be 
compared with the value cq + ^ log cr = 0.388(5). The estimates are compatible within the 
errors. Since the fit we performed depends on one free parameter only, our estimate for 
co + ^ log a turns out to be more precise than the one of [34] . As discussed in [34] the result 
we obtain is of the order of magnitude but somehow larger than the prediction obtained 
in [43, 39] in the framework of the (j)^ theory: cq + log cj ~ 0.29. 

To make the above discussion visually appreciable, in Figure 4 we plot our theoretical 
curve against the MC points, having fixed the value of the normalization M according to 
the fit at Lmin = 21 for the data set 1 and at Lmin = 11 for data set 2 and having truncated 
the sum over the oscillators at A^ = 100 in both cases. In fact, to draw a readable plot, it 
is convenient to subtract from both the Monte Carlo data and the theoretical prediction 
for the free energy the dominant scaling term a A and a constant^ term — ^ log a. 

In Table 9 of ref. [34] some data regarding rectangular lattices, i.e., with u > 1 were 
also presented. Such data represent an important test of the form of our expression 4.1. We 

This term is related to the natural appearance of a multiplicative factor a in the effective string 
partition function, see eq. (4.3). Furthermore we prefer to exhibit the same data plots appearing in [34]. 
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L2 


VaA 


u 


Fs 


difF {N = 100) 


difF {N = 0) 


10 


12 


2.29843 


6/5 


7.1670(6) 


0.0016 


0.0049 


10 


15 


2.56972 


3/2 


8.4449(12) 


-0.0004 


0.0042 


10 


18 


2.81498 


9/5 


9.6976(17) 


-0.0009 


0.0037 


10 


20 


2.96725 


2 


10.5235(25) 


-0.0012 


0.0035 


10 


22 


3.11208 


11/5 


11.3466(36) 


0.0017 


0.0064 



Table 2: Comparison between the data presented in [34] for rectangular lattices and our predictions. 
The data were extracted at /3 = 0.236025, corresponding to cr = 0.044023, i.e., in the same situation 
as for the square lattices considered in data set 2 considered above. In the last columns we give the 
MC values of the free energy Fs , according to Table 9 of [34] , and the difference between these values 
and the values obtained from our expression eq. (4.1), with the sum over the oscillator number 
truncated at A'' = 100 or at A'' = 0. The normalization constant J\f was already determined by the 
fit to the data set 2 presented in Table 1; we chose here the value obtained setting imin = 11- 



present the results of this test in Table 2; they are quite encouraging, since we basically get 
agreement with the data within the error bars, at least when including the contributions 
of higher oscillator numbers (up to N = 100). Notice that in this test we use the value 
of the normalization constant J\f already determined by fitting it to the data for squared 
lattices of data set 2, see Table 1, since the data for asymmetric lattices in ref. [34] were 
obtained at exactly the same value of the coupling constant f3 used for the data set 2. This 
makes the observed agreement even more remarkable, since it involved no adjustable free 
parameter. 

It would clearly be desirable to obtain more extended sets of data for asymmetric 
lattices, to check their degree of consistency with our Nambu-Goto predictions. 

One remarkable feature of the comparisons presented here in Tables 1 and 2 is that 
keeping just the zero modes of the string, i.e. setting N = 0, seems to yield a very good 
agreement with the data, in particular for symmetric interfaces. The contribution of non- 
zero modes modifies the overall normalization coming from the one-loop determinant, as 
discussed above, but it is much less important for higher loops, so that the shape of the 
free energy as a function of a A is little changed. 

Keeping only the zero-modes, the quantity £ of eq. (2.21) which appears in the parti- 
tion function reduces effectively to 



1 - "^(^-'^ (4.4) 
3aA ^ ' 



and the partition function eq. (4.3) simply becomes 



{^f\r^J'A' K._^{aAe). (4.5) 



5. Conclusions 

In this paper, we studied the free energy of interfaces in a compact d-dimensional target 
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space T"^, modeling them with a Nambu-Goto string effective action; this basicaUy corre- 
sponds, for the d = 3 case, to use the capillary wave model. Performing standard covariant 
quantization and explicitly integrating over the modular parameter of the world-sheet, we 
were able to derive an exact expression for the sector of the NG partition function that 
describes the interface fluctuations . 

Focusing on the three-dimensional case, we compared the exact predictions of the NG 
model with the free energy of interfaces obtained in [34] through very precise Monte Carlo 
simulations in the Ising model. The Monte Carlo data are very well accounted for when 
we consider lattices of sufficiently large sizes, typically with L > 2/-^^; here a is the string 
tension, which can be associated with great accuracy to the value of the Ising coupling [49] . 

Our analysis, which allows to compare the data with the exact theoretical prediction, 
confirms and strengthens the evidence that the NG model is a very reliable effective model 
for such lattice sizes. Previously, comparisons were made to the two-loop expansion of 
the NG model [29], which already gives a very good approximation for square interfaces; 
as pointed out in [34], however, and as confirmed here, higher order corrections are more 
important for asymmetric lattices. In fact, wc found extremely good agreement of the exact 
expression (which re-sums all the loop corrections) for the few data with u ^ 1 presented 
in [34]; it would be very interesting to obtain extended set of high precision Monte Carlo 
data for asymmetric lattices to check whether this agreement is confirmed, and to what 
extent. 

For smaller lattice sizes, rather large deviations from the NG predictions are observed 
in the data. This pattern is in perfect agreement with what has been observed in previous 
studies of the NG predictions for other observables, such as the Polyakov loop correlator 
[13], and is in fact to be expected for deep theoretical motivations. 

Indeed, we have studied the Nambu-Goto bosonic string in d dimension paying no 
attention to the fact that in a consistent quantum treatment, when d ^ 26 further degrees 
of freedom have to be taken into account, beside the d — 2 transverse oscillations. In the 
Polyakov first order formulation [46], the scale e'f' of the world-sheet metric h does not 
decouple, and gets in fact a Liouville-type action; in the Virasoro treatment, longitudinal 
modes of the string enter the game [50, 51] (the two effects are not disconnected, see [52]). 
For this reason these quantum bosonic string models have not been considered in the past 
[53, 57] as viable dual descriptions of situations, such as the confining regime of gauge 
theories or the fluid interfaces, where only transverse fluctuations are expected on physical 
grounds. 

In a modern perspective, the extra degree of freedom corresponding to the Liouville 
mode could play somehow the role of the renormalization scale [54, 55], in the spirit of (a 
non-conformal version of) AdS/CFT duality [56]. It would be very interesting to try to 
take into account the contributions to the partition function of the extra d.o.f. arising in 
d = 26 at the quantum level, and to see whether they can account for the deviations of 
data, at small lattice sizes, from the prediction obtained via the naive treatment of the NG 
model. Some effort in this direction, focusing on the relevance of the longitudinal modes 
for the description of large N gauge theories, has been done in [61]. 

In [57] Polchinski and Strominger put forward a different proposal to circumvent the 
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mismatch of d.o.f. between the bosonic string model and the transverse oscillations. An 
effective string was built, with an action analogous to the Polyakov one, where the inde- 
pendent scale e'l' of the metric gets replaced by the induced metric. This model is certainly 
very interesting, and it is receiving nowadays a renewed attention [58, 48]. It is noticeable 
that the excitation spectra obtained within this model agree with the naive NG ones up 
to two loops in the \/{aA) expansion; the prediction of this model would therefore nicely 
agree with NG in the region where the latter is valid. Higher corrections, whose compu- 
tation does not appear to be simple, could in principle better explain the data for smaller 
lattices. 

In conclusion, there is a renewed interest in the problem of identifying the correct QCD 
string dual and of testing the various proposals, sec for instance [19, 58, 59, 60, 48, 61]. 
As already said, our contribution corroborates the "experimental" result that the naive 
treatment of the NG model provides a very accurate description of fluctuating surfaces for 
sufficiently large sizes, and that clear deviations appear for smaller lattices which should 
be accounted for by the correct dual model (at least down to some lower scale where the 
string description itself breaks down). On the theoretical side, our derivation of the exact 
partition function, and not just of its perturbative expansion, could offer some insight for 
a similar treatment of consistent models including also the extra d.o.f. arising m. d ^ 26, 
or of the Polchinski-Strominger string. 
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A. Loop expansion of the exact result 

As discussed in section 2 eq.(2.23), the exact expression of the interface partition function 
x('^) = 2(^)~yTV^^^CfcC„-J^) ' K,_^{aA£) (A.l) 

m=Ofc=0 ^ ^ ^ 

has to be expanded for large aA, in order to compare it with the result of the functional 
integral approach. Therefore we have to use the asymptotic expression eq. (3.3) of the mod- 
ified Bessel functions for large argument which reads, making explicit the first coefficients 
af^ and ai\ 

^ y 2z \ 8z 2\{8zf 

Inserting the expansion of Bessel functions in eq. (A.l) gives eq. (3.4) (which we repeat 
here for commodity): 

I^'^ = V2^ {^)^ Vt Y: w(^) ' e-'^^[l + Y,ai'^\aA£)-A . (A.2) 

^ k,k'=0 V r=l / 
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The expression of £ was given in eq. (2.21), and can be rewritten as 



I 4tt u 47r2 ^2 



in terms of a and b defined as in eq. (3.6): 



d — 2 

a = k + k' — , b = k- k' . 



We expand £ as: 



£^a, 6) ~ 1 + —2na + j-j^2n' {b' - a') + j-j^^Tr^a {a' - b') + .. 



and insert this expansion into eq. (A. 2), obtaining 

-tA °° 



/ d-2 ' 
12 y 



x|l + -(d-2).a + ^(d-2) 



vr 



{d-6)—a^ + TrH'^ 



+ 



u 



l + ^2.2(a^-5^) 



ATT^u^a (a^ - 62) - 2n^u'' [b^ - a^f 
MP 



+ 



((7^)3 



X + 



1 d{d-2) 



d{d-2){rP-2(l-H) _ nv rH,d-2)a 
128 i 



+ 



1 



(tA 8 ((7^)2 

The result up to the third order is thus of the form of eq. (3.8): 



d-2 

^ ^ k,k'=0 



92 



aA {gAY 



with g\ and exphcitly given by 
d{d - 2) 



51 



92 



+ ^(d - 2)7ra - n227r2 {\?- - (?) 



d{d-2){d^ -2d- 
128 



Introducing the new variable Q = exp(— 2?™) we can write eq. (A. 4) in the form 



-a A 



T^'^ oc a——r^ lim <^ 1 



d-2 
U 2 Q^Q 



91 



+ 



92 



[riiQHQ)] 



2-d 



a A {aA)\ 

where we have replaced a and b with the derivatives on Q and Q (see eq. (3.10)): 



^ d d 
dQ dQ 
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As we can see from the expansion coefficients (A. 5) and (A. 6), the only terms involving 
powers of a and b higher than 1 are 



dQ dQ 

2 



dQ^dQ 



a^)' = 16 ( Q 



_d_ 

dQ 



dQ 



which have to be applied to rf '^{Q)'ff' '^{Q)- Applying the expressions eq. (C.8) and eq. 
(C.9) of the first logarithmic derivatives of the eta-function given in Appendix C, we have 



Q^n'-\Q) = {'^-dW-\Q)^ 



Q 



It is easy to check that, after the application of the derivatives, a factor of rf~'^'^{m) will 
always appear in front of the various terms and the total result up to the third order can 
be written in the following form: 



d-2 e 



1 



First of all we consider the /i term derived from the gi one. As one can see, the only 
derivative involved in the computation \s Q-^: this implies that our final formula will 
include E2 functions but no E^. We can divide the calculation in function of the powers 
of u. The term independent from it is the first one appearing in eq. (A.5): 

d{d - 2) 

/i,o - 

The next power of u, after the substitution (A. 8), gives, acting on the rj functions, 

h,i = -2iTu{d-2f 

The last term is instead: 



24 



/i,2 = 27rV(d-2) 



144 



The final result is: 
/i = 



(d-2) 



2 r 



(^\^El{iu)-'^uE2{:m 



+ 



d{d - 2) 



(A.IO) 



Proceeding in the same way to evaluate /2 we find: 



iT:\d-2f 
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(4 -d)El- 2E4 



-\ 2 



24 



+ u^^id - m - 2f 



^(4-d) 
12 ^ ^ 



E2E4 



+ u 



,7r\d-2f{d-4) 



E4 {3d-8)Ep 
72 144 



u^d{d - 2f{d - A)E2 + ^d{d - 2){d^ -2d -8) 



(A.ll) 
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(all modular forms above have to be evaluated at iu). 



B. The two loop contribution: functional integral computation 

In this appendix wc check the calculation of the two-loop terms of the free energy made by 
Dietz and Filk with the functional integral method. Starting from the partition function 
written in the physical gauge and expanding it in powers of one can evaluate it on a 
rectangular domain B of sizes {R, T) . Defining G as the inverse Laplacian of the theory, 
the second order terms are given, in the notation of [29], by: 



7(2) 



X,D{D - 1)((1) - (2)) + X2D{{S) + (4) + 2(1) - 4(2)) 



+ ^D(i^-l)((3) + (4)+6(l))] , (B.l) 

where Z? is d — 2 in our notation, F spcficics the topology of the boundary B; we are 
interested in the case where the topology is x Si. Ai and A2 are parameters which 
depend on the string model under consideration and in particular they are given by 

Ai = l, ^2 = -\ (B.2) 
in the Nambu-Goto case. Moreover, the following definitions are used: 

^'^ = ^Jo Jo ^'d^^'='''='dtdt' 

\z=z't=t' ) 



i2) = '-rdj\t(^ 

^Jo Jo \dzdt' 



(3) 



'0 ^0 
rR rT / 92Q ^ 2 



^' (TJo Jo \dtdt' 



T /52GX2 



\z=z',t=t' 



Notice that the formulsused by Dietz and Filk in eq. (3.4) have to be multiplied by — ^ 
to get the right normalization of their final result eq. (3.7). 

We shall deal with the divergent terms within the (^-function regularization scheme, 
following Dietz and Filk. In the Nambu-Goto model with x Si boundary conditions on 
the rectangular domain, the Green function is: 

1 exp{27ri^(t-t')}cxp{27ri|(z-z')} 

m,n= — oo ' 7^ 
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(2) 

Let us now evaluate the various terms of Zg^^g^ (rt)' begin with, we find 



(1) 



aRT 



E 



m.n,p.g— — oo 
(m.n)/(0.0) 
(p,9)/((),0) 



16 E y- 



+ (rn = 0, n / 0, j3, 



m,n,p,q=l \^-^ ^ W ) \W j 

+ (m, = 0, g 7^ 0) - (m = 0, n 7^ 0, p = 0, g 7^ 0)] 



(B.3) 



where with (m = 0,n ^ 0,p,q) we indicate the sum over any n different from zero, any p 
and any q when m has been fixed at zero. The last term in r.h.s. of eq. (B.3) is to avoid 
the double counting and the first one can be rewritten using the equality 



n' 



n 



m,n=l W ^ 



Eo i 



R\ TT R 



T 2AT 



(B.4) 



The result is 



+ {m,n,p = 0,q 0) - {m = 0,n ^ 0,p = 0,q ^ 0)] . 

The last three terms can be evaluated using (B.4) and the ^ function regularization 
and they sum up to 



---E, I i- 

6r 



T 



.R\ ttT 

1 — z:^ 

TJ 6R 

which vanishes using the modular property (C.6) for the Eisenstein series E2. We get thus 

1 



(1) 



-H 



aRT 

where we have defined for notational convenience 



Eo i 



R 



Eo i 



R 



(B.5) 



The contribution (2) vanishes because the terms in the sum are odd: 
(2) = 3^7; E -- = 



aRT ^ /"n^ , m2\ «2\ 

(m,n)#(0,0) V-" ^ J \"- 

(p,9)7^{0,0) 



To compute the expression ((3) + (4) + 6(1)) which appears in (B.l), one can rewrite 
it as a perfect square: 



((3) +(4) +6(1)) = - r dz r dt 

Jo Jo 



aRT 



dzdz' dtdt' 
2 



\z=z',t=t' 



ri 

m,n— — 00 R/ 
(m,n)5^(0,0) 



I lit 



+ 



T2" 



aRT 
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Collecting the previous results we obtain the following system of equations: 



aRT ' 
(2) = , 

(3) + (4) + 6(l) = 
which gives, for the last term in equation (B.l): 

O ATT 

(3) + (4) + 2(1) - 4(2) = . 
The total second order correction is therefore 

Zs,x»„B = -^iD(D-l)H-^D{3-iH)--^D{D-l) 

This corresponds to our result^ eq. (A. 10) where D = d — 2, up to an overall factor of —5. 



C. Useful formulae 

In this appendix we collect some useful formulae. 

Dedekind rj function The Dedekind eta function is defined, in terms of the quantity 
q = exp{27rir}, by 

00 

V{r) = n (1 - ■ (C.l) 

n=l 

One can expand it in g-series: 

ivir)]-' = ^Pkq'-^^ , (C.2) 
fc=o 

where p^.'s arc the number of partitions of k. In the text, wc often switch between the 
notation r]{T) and r]{q) for the function defined in cq.s (C.1,C.2), according to the conve- 
nience. 

Under the modular transformations T and S the Dedekind eta function transforms in 
the following way: 

ry(r + l) = eitry(r) , (C.3) 
77(--) = {-iT)-2r,iT) . (C.4) 



*Dietz and Filk found: zf^^^ - [D^H - i[/)(4D - 1)]} 
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Eisenstein series The second Eisenstein function is defined by 



E2{t) = 1 - 24^ai(n)g" = 1 - 24^ ^ 



n=l k=l ^ 

where ai{n) denotes the sum of the positive divisors of n. An useful property is 

where the fourth Eisenstein series is defined, in terms of the sum of the cubes of the positive 
divisors of n, (T3(n), as 

+ 00 

E4{r) = 1 + 240^(73 (n)g" . 

n=l 

The modular properties of the Eisenstein functions are: 

^ / 1\ 6 i 



E2(.r)=( -) E,^--j+--, (C.6) 



£^2fc(T) = (-lf - £;2fe— ,fe>2. (C.7) 



Multiple logarithmic derivatives of the eta function are related to Eisenstein series, and in 
particular we have: 

,^,"(T) = arf(r)^ , (C.8) 

Poisson resummation formula In section 2 the following resummation formula plays 
a key role: 

+ CXD 

I 'TTi m — n I- I 

(C.IO) 



Er 2 « 1 -i I Trim — 0) \ 
expj— Tran + 2mbn\ = a 2 \ exp < > 
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